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In this paper we study how to attack under the self-similarity hypothesis a perfect fluid 
Bianchi I model with variable G,and A, but under the condition divT ^ 0. We arrive to 
the conclusion that; G and A are decreasing time functions (the sing of A depends on 
the equation of state) , while the exponents of the scale factor must satisfy the conditions 
X]i=i o^i = 1 3.nd '^i < li G 1) 1 relaxing in this way the Kasner conditions. 

We also show the connection between the behavior of G and the Weyl tensor. 

Keywords: Time varying constants; Bianchi I; Self-similarity. 



1. Introduction. 

Ever since Dirac first considered the possibility of a G variable (see ^ , there have 
been numerous modifications of general relativity to allow for a variable G, nev- 
ertheless these theories have not gained wide acceptance. However, recently (see 
1^11131) a modification has been proposed treating G and A as non-constants coupling 
scalars. So it is considered G and A as coupling scalars within the Einstein equa- 
tions, Rij — — GTij — A^y , while the other symbols have their usual meaning 
and hence the principle of equivalence then demands that only gij and not G and 
A must enter the equation of motion of particles and photons. In this way the usual 
conservation law, divT = 0, holds. Taking the divergence of the Einstein equations 
and using the Bianchi identities we obtain the an equation that controls the vari- 
ation of G and A. These are the modified field equations that allow to take into 
account a variable G and A. Nevertheless this approach has some drawbacks, for ex- 
ample, it cannot derived from a Hamiltonian, although there are several advantages 
in the approach. 

There are many publications devoted to study the variation of G and A in the 
framework of flat FRW symmetries (see for example ^^IS-j g^j-j^j ^}jjg -^yorks have 
been extended to more complicated geometries, like for example Bianchi I models, 
which represent the simplest generalization of the flat FRW models (see for example 
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1^^^ . in the context of perfect fluids and ^^^^ in the context of viscous fluids) . 
Bianchi I models are important in the study of anisotropies. 

But in our opinion, the problem arises when one try to solve the resulting field 
equations (FE). It seems that it is unavoidable to make simplifying hypotheses, or 
to impose ad hoc some particular behavior for some of the quantities of the model, 
in order to obtain a exact solution to the FE. Such simplifying hypothesis are made 
for mathematical reason (in order to reduce the number of unknowns) although are 
justified form the physical point of view. Usually such assumptions or simplifying 
hypotheses follow a power law, for example, the quantity X follows a power law i.e. 
X = Xot" , where Xq is an appropriate dimensional constant, t is the cosmic time 
(for example) and a G M (usually a G Q, but this other question), and depending 
on the nature of the quantity X, a will be positive or negative. Actually we think 
that although all these simplifying hypotheses are correct or at least bring us to 
obtain correct results, it is not necessary to do that, since they may be deduced 
from symmetry principles in such a way that one may justify (deduce) them from 
a correct mathematical principle, and usually all these approaches have physical 
meaning. 

Therefore the main goal of this paper is to apply the well known tactic (approach) 
of self similarity (SS) in order to study and find exact solutions for a perfect fluid 
Bianchi 1 models with variable G and A, but under the condition divT ^ 0, and 
trying to make the lowest number of assumptions or neither. We will try to show 
that with this approach all the usual simplifying hypotheses may be deduced from 
a correct mathematical principle. 

The paper is divided in the following sections. Section two is devoted to outline 
all the ingredients as well as the field equations. In section three, we study the 
model under the self-similarity hypothesis. We start this section introducing briefly 
some ideas about self-similarity and self-similar spacetimes. Once we have found 
the homothetic vector field we go next to calculate the scale factors, where they 
obviously follow a power law solution, as well as the derived quantities from them 
as the Hubble parameter H, the deceleration parameter q and the shear ct, since 
they only depend of the scale factors and will be the same for all the studied cases. 
We study four models. The first of them is the classical one i.e. where G is a 
true constant and the cosmological constant A vanish. We have preferred starting 
with this model in other to check how works the employed tactic. We emphasize 
that the obtained solution (that will be the same in all the studied cases) satisfies 
the condition, X]i=i ~ 1 ^^'^ X]i=i o^l 1' relaxing in this way the Kasner 
conditions {J2i=i'^i — J2i=i'^i — 1j only valid for a vacuum solution), but it is 
only valid for uj — 1, i.e. ultra-stiff matter, and where (ai)^=i ^"^^ the exponents of 
the scale factors. liuj ^ 1, the model collapse to the standard flat FRW one. We then 
study the curvature invariants as well as the Weyl tensor and its invariant, and end 
calculating the gravitational entropy. We show that the definitions of gravitational 
entropy do not work well in this kind of spacetimes, the self-similar ones, since these 
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definitions are dimensionless which means that this quantity remains constant along 
the homothetic trajectories of any self-similar spacetime. We would like to relate in 
same way the behavior of the Weyl tensor with the behavior of a G time- varying. 
So after showing that the tactic works correctly, we pass to study a simple model 
where we consider that only vary G. In this case we get the same solution as in the 
above model with the same restriction for the equation of state i.e. the solution is 
only valid ii lj = 1. In the third of the studied model we consider only a A time- 
varying. As we have pointed out previously the solution is the same with regard to 
the exponents of the sale factors i.e. they must satisfy the conditions X]i=i c^i = 1 
and Y^i=i '^i < 1' t)ut in this case this solution is only valid iff w G (—1, 1) i.e. uj 1, 
since ii ui — 1, A — 0. We find always that A behaves as a negative decreasing time 
function. The last of the studied cases considers that both "constants" , G and A 
are time-varying. We find that G is a decreasing time function (i.e. has a similar 
behavior as the Weyl tensor) and A « t^^, with A < 0, and as in the above case 
this solution is only valid Vw G (—1,1). li uj > 1 then A > 0. We also show how to 
regain the "classical" case where div T = 0. We end we a brief conclusions. 

2. The Model. 

Throughout the paper M will denote the usual smooth (connected, Hausdorff, 
4-dimensional) spacetime manifold with smooth Lorentz metric g of signature 
(— , -I-). Thus M is paracompact. A comma, semi-colon and the symbol C 
denote the usual partial, covariant and Lie derivative, respectively, the covariant 
derivative being with respect to the Levi-Civita connection on M derived from g. 
The associated Ricci and stress-energy tensors will be denoted in component form 
by Rij{= R'^jcd) and Tij respectively. A diagonal Bianchi I space-time is a spatially 
homogeneous space-time which admits an abelian group of isometrics G3, acting on 
spacelike hypersurfaces, generated by the spacelike KVs ^1 = dx,S,2 = dy,£,3 = dz- 
In synchronous co-ordinates the metric is: 



where the metric functions Ai{t), A2(t), A3{t) are functions of the time co-ordinate 
only (Greek indices take the space values 1,2,3 and Latin indices the space-time 
values 0,1,2,3). In this paper we are interested only in proper diagonal Bianchi 
I space-times (which in the following will be referred for convenience simply as 
Bianchi I space-times), hence all metric functions are assumed to be different and 
the dimension of the group of isometrics acting on the spacelike hypersurfaces is 
three. Therefore we consider the Bianchi type I metric as 



ds^ = -dt^ + Al{t){dx^')^ 



(1) 



c^dt^ + X^{t)dx^ + Y^{t)dy^ + Z^{t)dz^, 



(2) 




For a perfect fluid with energy-momentum tensor: 



Tij ^ (p + p) u^Uj + pgij , 



(3) 
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where we are assuming an equation of state p — top, [oj = const.). Note that here 
we have preferred to assume this equation of state but as we will show in the 
following sections this equation may be deduced from the symmetries principles as 
for example the self-similar one. The 4— velocity is defined as follows 

u = (^i, 0,0,0^ , um^ = -l. (4) 

The time derivatives of G and A are related by the Bianchi identities 

1„ /87rG' ^'^ 

2 



R^3 - o^ff^j = —T^iJ - ^9i] ) , (5) 



in this case this equation reads: 

Therefore the resuhing field equations are: 

—p + Ac", (7 
cr 

-uop + Ac^, (8 
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ujp + Ac^, (9) 



cjp + Ac^ (10) 

-8^-G^- 
Now, we define 

Since we have defined the 4-velocity by eq. ^ then the expansion is defined 
as follows: 

» = l(|4 + |).iH, ,13) 

and therefore the acceleration is: a.^ = Ui.jU^,m this case a = 0, while the shear is 
defined as follows: o-y = ^ {ui-ahj + Uj-^ah^) — ^^hij, 




2 1 9 1 / / a: \ / y \ / z \ XY xz 
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3. Self-similar solution. 

In general relativity, the term self-similarity can be used in two ways. One is for the 
properties of spacetimes, the other is for the properties of matter fields. These are 
not equivalent in general. The self -sim ilarity in general relati vity wa s defined for the 
first time by Cahill and Taub (see^SS^ and for general reviews ^^SUlll). Self-similarity 
is defined by the existence of a homothetic vector V in the spacetime, which satisfies 

£v5y = 2a5ij, (15) 

where gij is the metric tensor, Cv denotes Lie differentiation along V and a is a 
constant. This is a special type of conformal Killing vectors. This self-similarity is 
called homothety. If a ^ 0, then it can be set to be unity by a constant rescaling of 
y. If a = 0, i.e. Cygij = 0, then is a Killing vector. 

Homothety is a purely geometric property of spacetime so that the physical 
quantity does not necessarily exhibit self-similarity such as CyZ = dZ, where d is 
a constant and Z is, for example, the pressure, the energy density and so on. From 
equation it follows that 

CvR'jki=0, (16) 

and hence 

CvR^J = 0, CvG^j = 0. (17) 

A vector field V that satisfies the above equations is called a curvature coUineation, 
a Ricci coUineation and a matter collineation, respectively. It is noted that such 
equations do not necessarily mean that F is a homothetic vector. We consider the 
Einstein equations 

Gy=8^GTy, (18) 

where Ty is the energy-momentum tensor. 

If the spacetime is homothetic, the energy-momentum tensor of the matter fields 
must satisfy 

CvT^j = 0, (19) 

through equations p8)l and (|17p . For a perfect fluid case, the energy-momentum 
tensor takes the form of eq. ([3]) i.e. Ty = {p + p)uiUj -|- pg^j , where p and p are the 
pressure and the energy density, respectively. Then, equations p5|) and ^9]) result 
in 

Cyu^ = —au^, Cvp = ~'2ap, Cyp = —2ap. (20) 

As shown above, for a perfect fluid, the self-similarity of the spacetime and that 
of the physical quantity coincide. However, this fact does not necessarily hold for 
more general matter fields. Thus the self-similar variables can be determined from 
dimensional considerations in the case of homothety. Therefore, we can conclude 
homothety as the general relativistic analogue of complete similarity. 
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From the constraints (|20p , we can show that if we consider the barotropic equa- 
tion of state, i.e., p = f{p), then the equation of state must have the form p — up, 
where is a constant. This class of equations of state contains a stiff fluid {uj = V) 
as special cases, whiting this theoretical framework. There are many papers de- 
voted to study Bianchi I mode ls ("in d ifferent context) assuming the hypothesis of 
self-similarity (see for example SSMZl) here, we would like to try to show how 
taking into account this class of hypothesis one is able to find exact solutions to the 
field equations within the framework of the time varying constants. 

The homothetic equations are given by eq. (|15p so it is a straightforward task 
to find the homothetic vector field, where in this case is as follows: 

V^tdt+{l- xd, + (^1- t^^ ydy zd,, (21) 

Therefore, we have obtained the following behavior for the scale factors: 

X = Xo^i , Y = Yot"' , Z = Zot""' , (22) 
with Xq, Yq, Zq are integrating constants and (ai)^^-^ G M. In this way we find that 





(23) 

with (ai + a2 + as) = a. 

In this section we are going to study several Bianchi I models and we will show 
how it is possible to find exact solutions to the field equations (without the condition 
divT = 0) under the hypothesis of SS. 

The time derivatives of G and A are related by the Bianchi identities i.e. eq. 
([TT|) that in this case collapses to the following one: 

p + p{l + c.)H = fit) = - -p, (24) 

where f{t) is a function that depends on time and controls the time variation of 
the constant G or/and A. If G = const, and A vanish then f{t) = 0, so the model 
collapses to the standard one. This idea was pointed out by Rastal (see 1^ and 
improved (in the theoretical framework of time varying constants) by Harko and 
Mak (seel^^. 

Therefore the resulting field equations are ([7lfT0|) together to the new one 

p + p{l+uj)H^f{t), (25) 

3.1. "Constants" constants. The classical model. 

In this case we consider f(t) = 0, so this means that G — const, and A vanish and 
therefore we get that from eq. ((25|) that 

p + pil+Lj)H ^0, p = poi"^"+'^"- (26) 
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From the field equations ([7]) we get that 

Ac^ 2 

where A = aia2 + a^ai + 0203. 

The shear has the following behavior, cr^ ^ 0, as it is observed cr — > as 
{ai —>■ aj) . As in the previous sections, we may calculate the coefficients (a^) by 
solving the following system of equations: 

"2 (a2 - 1) + "3 (0:3 - 1) + Q:3Q!2 = -Auj, (28) 

ai (cti ^ 1) + 0:3 (0:3 ^ 1) + ct3Q!i = —Auj, (29) 

a2 (a2 — 1) + ai (ai - 1) + 0102 = -At^^, (30) 

a{uj + l) = 2, (31) 

where A — a\ai + asai + 0203, and a = ai + 02 + ck3- 

So we have the following solutions for this system of equations: 

ai = l-a2-a3, = 1, (32) 

ai=a2 = Q;3 = — — — , (33) 

0\L0 ^ I) 



as it is observed only solution ([32)1 is interesting for us. The second solution is the 
usual FRW one, so it is not interesting for us (see Einstein&de Sitter (j^^ for w = 0, 
and Harrison ^'^^ Vw). Nevertheless we have found that solution ([5^ verifies the 
conditions 

a = = 1, ^a\<\, (34) 

but iff w = 1 , (see l^^h^Jiiie other authors claim that must be satisfies the con- 
dition 'Yli'^i = 1; (see SH Bni and and in particular, in this context (seeS^. 
Nevertheless we have found that this solution only verifies the first of the condi- 
tion of the Kasner like solutions i.e. a = a,- = 1, while the second condition 
^ a| — 1, it is not verified (see ^^Hl and^^Sl), this case we find that it is 
verified the conditio n V] af < 1. Therefore we have found the same behavior as 
the obtained one in p^. Before ending we would like to make a little comment 
about the Kasner like solutions. If a solution of (|28ll3ip verifies the relationships 
2i Q!? = X^^tti = 1, i.e. they are Kasner's type (see 42,^ 31! and in particular , 
then this means that A = a\OLi + olj.olx + ol20lz — (i.e. the model is Ricci flat), 
which brings us to get the following result: ai = i ^1 — Q!3 — ^1 + 2a3 — 3a§^ < 0, 

a2 = 5 1^1 — 0:3 + v^l 4- 2a3 — > 0, Va3 G (0, 1) , we think that this class of 

solutions are unphysical and have a pathological curvature behavior as it is shown 
bellow. Furthermore, as we can see, if A = 0, then from eq. ([27|l we get p = 0, as it 
is expected for this class of solutions (vacuum solutions) so they are not interested 
for us. Nevertheless relaxing the condition J^'^f — to our result i.e. X)*^? < 1j 
we are able to obtain solutions whit (a.;) > 0,Vi, and p ^ 0. 
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Therefore we have obtained the following behavior for the main quantities: 

H^^, => g = 0, (35) 

so it is quite difficult to reconcile this model with the observational data. With 
regard to the energy density we find that 

P-^t-^ -^-^(1 + 3A)1, (36) 

and with regard to the constants (cui)'^^^ we have that only obtain a BI solution 
iff ai — \ — a2 ~ as, (where furthermore we suppose that ck2 ^ Q^s) and that this 
res ult only is possible if the equation of state is w = 1, i.e . ultra-stiff matter (see 
ESI) . For a review of Bianchi I solutions see for example (j^^ . 

With regard to the curvature behavior, we may see that the full contraction of the 
Ricmann tensor (see for example 1^1^^ Ii ■= Rijj^iR^^^^ ^ and the full contraction 
of the Ricci tensor, I2 := RijR^\ are: 

K 'i{-a2 - as. + al + a2az + alf 
h - h- , (37) 

where K = K{ai) = const ^ 0, i.e. 

K =[i {al + al) + 2a2a3 + 9a^a^ + 3 (a^ + af) - 6 (a^ + al) + 02^3 (6 {al + a^) - 8 (02 + 03))] . 

(38) 

and the scalar curvature i? is: R — R\. 

The non-zero components of the Weyl tensor are: 

C1212 = K^t-^^^^^'\ Cisrs = K2t-^^^-'''\ Ci4i4 - i^3t"'('~"^\ 

C2323 = ^4^~'"^ ^2424 = C3434 = i^6i"'"\ (39) 

where the numerical constants (ifi)i'=i — K{ai) — const ^ 0. As we can see with 
the obtained solution for {oLi)^^-^ the Weyl ^ 00 as t — > 0, in the next models we 
shall show that G{t) has the same behavior as the Weyl tensor. In a forthcoming 
paper we study models with Weyl ^ as t — > (i.e. models that verify the Weyl 
tensor hypothesis) and with a growing G, i.e. in same way exists a relationship 
between both quantities. 

Now taking into account a very famous result by Hall et al (see S3) we may 
check that CyC^ju = 0, as it is shown in (S^ if a vector field V G X(M), verifies the 
conditions CyCji^i — 0, and CyTij = (as it is known, if y is a homothetic vector 
field, then it is also a matter coUineation), then Cyg — '2g i.e. it is a homothetic 
vector field, but in this case we have arrived to the conclusion that Cyg = 2^ 
CyTij — 0, and that it is also verified the relationship LyCji^i — 0. 

The Weyl scalar is defined as: I3 :— C"''"^'^Cabcd, as it is observed /3, is also 
defined as follows: I3 — Ii — 2/2 -I- j-R^, this definition is only valid when n = 4. 
Therefore, I3 has the following behavior 

^3 = (40) 
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with K given by K = 

^ [al + ^3 - a2a3 + -iajal + aj + al-2[al + a^) + a2a3 (2 (a| + af) - {a2 + as))] 
The non-zero components of the electric part of the Weyl tensor are: 

^22 - ^33 = if2<"'^'~"'\ ^44-^^3^2(1-"^), (41) 

while the magnetic part of the Weyl tensor vanish, Hj j = 0. 
The gravitational entropy is defined as follows (see I '^^H'^^ I) : 

„ /g h- 2/2 - I ir2 

P' = ^ = - , ' =7^+o^~2. (42) 



h h h 3 1: 



2 



finding that 



P'' = const. ^ 0, (43) 

note that = h/h- As have been pointed out by Nicos Pelavas et al (see^^^l) this 
definition is not an acceptable candidate for gravitational entropy along the homo- 
thetic trajectories of any self-similar spacetime. Nor indeed is any "dimensionless" 
scalar. This implies that I3/I2 is constant along timelike homothetic trajectories. 
As a consequence, (j42p does not provide a measure of gravitational entropy along 
homotheticities and therefore /3//2 cannot be a candidate for a measure of gravi- 
tational entropy in self-similar spacetimes. 



3.2. G— variable. 

In this case we are going to consider that only vary "constant" G. This only possible 

G ' 



if we take into account the condition divT ^ and therefore f{t) — — so eq. 



([25)1 collapses to the following one. 

^ + ^ = -(1+^)7, =^ pG^t-('+-^", (44) 
p G t 

From the field equations ([7]) we get that 

A 1 2 
^P^i^^lT^¥^ ^^^^ 

The shear has the following behavior, ^ 0, as it is observed cr ^ as 
(tti aj) . As in the previous sections, to calculate the coefficients (a^) we need to 
solve the filed equations obtaining the same system of equations as in the above case 
i.e. eqs. (|28ll3ip . so we get the same solution as in the above case i.e. X^i'^i '^i ~ 1' 
and ^1 1' ^ii'^ ^ii^y valid if w = 1. 

Therefore we have obtained the following behavior for the main quantities: 

H^\. =^ g = 0, (46) 

and 
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note that this result is quite similar to the obtained one in the last solution i.e. the 
obtained one in eq. ((36)) . but we are not able to get a separate behavior for the 
quantities G and p. 

3.3. A— variable. 

In this case we consider only the variation of the cosmological constant A, so eq. 
((25)) yields 

and therefore from the field equations ((7)) we get that 



and hence 

c2 A 1 



^ 47rG (1 +w)a^2■ 
Now, we next to calculate the quantity A, from eq. ((49)) we get 



(50) 



A^4fl- .. \ (51) 
c2 V (1 + w) a/ i2 ' V ^ 

in this way it is observed that 

r Ao > <^ {uj + l)a>2 
A = Ao^"^ < Ao = ^ (w + l)a = 2 . (52) 
i Ao < {uj + l)a < 2 

The shear has the following behavior, ^ 0, by hypothesis. As in the previous 
sections, we calculate the coefficients (cti) from following system of equations: 



0L2 {oL2 - 1) + as (as - 1) + a3a2 = A \^ ^ j > (53) 
ai (ai - 1) + as (as - 1) + asai = A ^ i (54) 
Oi-2 [oL-i - 1) + ai (ai - 1) + aia2 = A \- — ^ ) , (55) 



where A — aia2 + asai + a2a3, and a — ai + a2 + a^. 

So we have the following solutions for this system of equations: 

ai = a2 = as, (56) 
ai = 1 - a2 - as, (57) 

as it is observed solution (|56p is not interesting for us, since it is unphysical (in this 
context). Only the second solution has physical meaning and it is valid G ( — 1,1]. 
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Therefore, we have found that the solution is Y^i=i '^i ~ 1' Si=i '^i 1' 
in occasion this solution is only valid if Vw g (—1, 1] . 

Therefore we have obtained the following behavior for the main quantities: 

H^j, ^ g-0, (58) 



7^-T72' (59) 



while with regard to the energy density we find that 

A 1 

AttG {1 + lo)T^ 

so, if ti> < — 1 =^ p is negative (phantom cosmologies), for the rest of the values of 
LU, i.e. to G (—1, 1], p is a decreasing function on time. 
The cosmological "constant" behaves as follows 

^-^ot-^ ^o = \^;=l^^=l\ , , (60) 
1^ Ao < 0, V e (-1, 1) 

so we have found that A is a ^^negative decreasing function^^ on time. Note that 
Ao > iff w > 1. As we can see this solution is quite different of the previous ones, 
since here we have obtained a solution type Bianchi I V uj £ (—1,1) while in the 
previous ones this only happens if lu = 1. Here if cj = 1 then we regain the first of 
the studied cases i.e. which one where A vanish and G behaves as a true constant. 

3.4. G&iA— variable. 

In this case we are going to consider that both "constants" G and A vary, therefore 
eq. ([25]) yields 

p + p(l+cu)H = -^-^p, (61) 
ottG G 

so from the field equations (O and (fTTj) and following the same steps as in the 
above models we get that 

Gp=f^^—t-\ (62) 
47r (UJ + l)a 

as we can see it is verified the relationship Gp « t~'^, as it is expected. In fact 
it is impossible to separate both functions (to get the behavior of both functions 
independently), to do that we need to impose a condition, but precisely we are 
trying to avoid such way. 

Now taking into account again eq. ([7]) we get 

A = Ao^-^ Ao^^fl--^—) , (63) 

\ [LU + l)a J 

in this way it is observed that 

{lu + l)a > 2 

A = Ao^"^ { Ao = Q ■iF=^ {uj + l)a = 2 . (64) 

{uj + l)a < 2 
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The shear behaves (see eq. (fT4|)) as foUows: ^ 0, by hypothesis. In order to 
find the value of constants (0;^) , we make that them verify the field eqs. so in this 
case we get the following system of eqs.: 



a2 {a2 - 1) + as [as - 1) + "3^2 = A y ^ j 1 (65) 
ai (ai - 1) + as (03 - 1) + as"! = A (^—^^ > (66) 
a2 {a2 - 1) + ai (ai - 1) + aia2 = A (- — -] , (67) 



(note that this system is the same as eqs. ()53ll55p where A = aia2 + asai + a2a3, 
and a = Q!i + a2 + ck3. Therefore we obtain the same solution as in the last studied 
case, i.e. Q^i = li ^-nd "^f < li with w £ ( — 1,1], finding in this way the 

following behavior for the main quantities: 

= q^O, (68) 

and with regard to the product Gp we get 

^^-£(^^"^' 

but we cannot say anything more. The cosmological "constant" behaves as follows 

so we have found that A is a negative decreasing function on time. As in the above 
case we get a positive cosmological "constant" if w > 1. 

In order to try to find a separate behavior for the functions p and G, we may 
suppose that 

P = p,t--, G = Got^-', =^ Gp^^-^t-^ =Kt-', (71) 

in [LU + 1) 

with a G K^, i.e. for example we may choice 

r'^ A K 

G = / , = -i""', (72) 

47rpo(w + l) Po 

therefore, it is verified the field eq. ([6T|) for all the possible values of a. We may find 
other possibilities as for example 

P = Po<"'+^ G = Go<-'^ (73) 

with a G (—00,2) . At this point we would like to stress the relationship between 
the behavior of the Weyl tensor and the behavior of G(i), since Weyl — > 00 as 
t ^ in the same way as G{t). 
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But if a = (ti> + 1) then we regain the condition divT = as well as f{t) — 0, 

diyT^p + p{l + iu)H = = -—--p = f{t). (74) 
but this case has been already studied in (j^^. 

4. Conclusions. 

We have shown how to attack a perfect fluid Bianchi I with G and A variable under 
the condition divT ^ 0, taking account only the hypothesis of SS. Our arguments 
exploit the symmetry properties of homothetically self-similar spacetimes. These 
calculations are of physical interest since self-similar spacetimes are very widely 
studied ^21 and are, for example, believed to play an important role in describing 
the asymptotic properties of more general models. In this way we have shown that 
it is not necessary to make any assumption "ad hoc" or to take into account any 
previous hypothesis or considering any hypothetical behavior for any quantity since 
all these hypotheses could be deduced from the symmetry principles, as for example 
the self-similar hypothesis SS (for other approaches see^i%. As have seen to get the 
solution under the SS hypothesis it very simple since once we have calculate the 
homothetic vector field it is a trivial task to obtain the behavior of the scale factors, 
where the must follow a power law solution, X = XqI"^, Y = Yot"^ , Y — YqI"^ , in 
this way we only need to calculate the value of theses exponents i.e. (ai)j^]^ , since 
the rest of the quantities will depend on these values. We have started studying the 
simplest case where G is a true constant and A vanish, i.e. the classical model in 
order to check how works the purposed method. For this model we arrive to the 
already known result = 1, and < 1, stressing that this result is only 

valid if the equation of state verifies the relationship w = 1, i.e. the result is only 
valid for ultra-stiff matter, otherwise the model collapses to the FRW solution. We 
have discussed why it is not possible to get the vacuum {p = 0) Kasner solution, 
^ai = l = ^Q!f = l. We think that this class of solutions are unphysical since 
necessarily one of the scale factors must be a decreasing time function (maybe such 
class of solutions would have any interest in the study of singularities). At the same 
time, we have studied the curvature invariants, Ii and I2, i.e. the Kretschmann's 
scalars, showing that the obtained solution is singular as well as the Weyl tensor 
and its scalar, I3. We have performed all these calculations in order to show that 
there is a relationship between the behavior of the Weyl tensor and the behavior of 
the variable Newton constant G{t). All these considerations are valid for the rest of 
the studied models. We also have calculated the gravitational entropy, — I3/I2, 
showing that this definition is not valid for self-similar spacetimes since this quantity 
is dimensionless and therefore remains constant. Actually we h ave arrived to the 
same conclusion as the obtained one by Nicos Pelavas et al (see^^^. Furthermore 
if one gets the Kasner solution then I2 ~ 0, i.e. the model collapses to a Ricci flat 
model and therefore = 00. 
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With regard to the second of the studied models, where only vary G{t), we have 
shown that it is not possible to get a separate behavior of the "constants" G and 
/5, obtaining Gp t^^, as it is expected. For this model we have obtained the same 
solution as in the previous case i.e. that the exponents of the scale factors must 
satisfy the relationships = 1, and ^af < 1, and only valid for w = 1. As we 

have shown in the last of the studied models in order to get a separate behavior 
between G and p, it is possible to follow several ways but in this paper we are trying 
to do the slightest number of hypotheses and in any case to avoid to make previous 
assumptions on the behavior of any quantity. 

In the third of the studied models, where only vary the cosmological constant A, 
we have found that A is a negative decreasing function on time Vw G (—1,1) and a 
positive time decreasing function if a; > 1. If w = 1 then A vanish, so the exponents 
verify the same relationship as above but in this case this solution is only valid if 
u! G (—1,1) since ii uj = 1 then we get the first of the studied models. 

In the fourth model, we have considered that both constant vary. In this case 
we arrive to similar conclusions as in the above cases, i.e. A is a negative decreasing 
function on time \/lu G (—1, 1), vanish if = 1 and it is a positive decreasing time 
function if > 1 while Gp ^ t^"^ and the exponents of the scale factors must 
satisfy the relationships "^ai = 1, and Y^af < 1, valid G (—1,1). In this 
occasion we have made an assumption on the behavior of G in of to try to know 
its behavior finding in this way that is a decreasing time function on time. It is 
quite surprising result since in similar models with FRW symmetries this quantity 
is always growing. At the same time we have shown the similitude between its 
behavior and the behavior of the Weyl tensor. Both quantities tend to infinite as t 
runs to zero We have finished howing how to regain, in a trivial way, the condition 
divT = 0. 
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